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^ ■ Abstract. 

^ ' We present a review on the determination of the primordial helium abundance Yp, based on 
the study of hydrogen and helium recombination lines in extragalactic H ll regions. We also 
discuss the observational determinations of the increase of helium to the increase of oxygen 
by mass AY/ AO, and compare them with predictions based on models of galactic chemical 
, evolution. 
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Q 1. Overview 

, The determinations of the hehum abundance by mass, Y, from metal-poor extragalac- 
^ ' tic H II regions provide the best method to obtain the primordial helium abundance. 
I— I , During their evolution galaxies produce a certain amount of helium and oxygen per unit 
mass that we will call AY and AO. The galaxies less affected by chemical evolution are 
^ , those that present a large fraction of their baryonic mass in gaseous form and a small 
' fraction of their baryonic mass in stellar form. These galaxies when experiencing bursts 
\^ of star formation present bright metal poor H 11 regions that have been used to determine 
their chemical composition. 

From a set of Y and O values and assuming a linear relationship it is possible to obtain 



o 



Yp and AY/ AO from the following equation: 



.AY 



Yp-Y-0—. (1.1) 



^ ■ The determinations of 1^ and AY/ AO are important for at least the following reasons: 
k>( (a) Yp is one of the pillars of Big Bang cosmology and an accurate determination of 
j_j Yp permits to test the Standard Big Bang Nucleosynthesis (SBBN), (b) the models 
^ of stellar evolution require an accurate initial Y value; this is given by Yp plus the 
additional AY produced by galactic chemical evolution, which can be estimated based 
on the observationally determined AY/ AO ratio, (c) the combination of Yp and AY/ AO 
is needed to test models of galactic chemical evolution. 

Recent reviews on primordial nucleosynthesis have been presented by Steigman (2007), 
Olive (2008), Weinberg (2008), and Pagel (2009). A review on the primordial helium 
abundance has been presented by Peimbert (2008) and a historical note on the primordial 
helium abundance has been presented by Peimbert & Torres-Peimbert (1999). 



1 



2 



Peimbert et al. 



2. Recent Yp determinations 

The best Yp determinations are those by Izotov, Thuan, & Stasinska (2007) that 
amounts to 0.2516 ± 0.0011 and Peimbert, Luridiana, & Peimbert (2007a) that amounts 
to 0.2477 ± 0.0029. The procedures used by both groups to determine Yp are very dif- 
ferent and it is not easy to make a detailed comparison of all the steps carried out by 
each of them. We consider that the error presented by Izotov et al. is a lower limit to the 
total error because it does not include estimates of some systematic errors. The difference 
between the central Yp values derived by both groups is mainly due to the treatment of 
the temperature structure of the H II regions. Izotov et al. adopt temperature variations 
that are smaller than those derived by Peimbert et al. 

To be more specific we can define the temperature structure of the H ii regions by 
means of an average temperature, Tq, and a mean square temperature fluctuation, 
(Peimbert 1967). The value of derived by Izotov et al. (2007) for their sample is about 
0.01; while Peimbert et al. (2007a) obtain a t'^ of about 0.026. From the observations 
adopted by Peimbert et al. and assuming = 0.000 we obtain Yp = 0.2523 ± 0.0027, 
and for t"^ = 0.01 we obtain Yp = 0.2505. 

The small value of t"^ derived by Izotov et al. (2007) is due to the parameter space 
used in their Monte Carlo computation where they permitted T(He i) to vary from 0.95 
to 1.0 times the T(4363/5007) value, which yields a of about 0.01. By allowing their 
r(He II) to vary from 0.80 to 1.0 times the T(4363/5007) value their t'^ result would 
have become higher. This can be seen from their Table 5 where 71 of the 93 spectra 
corresponded to the lowest T(He i) allowed by the permitted parameter space of the 
Monte Carlo computation. A higher t"^ value for the Izotov et al. (2007) sample produces 
a lower Yp, reducing the difference with the Peimbert et al. (2007a) Yp value. 

Two other systematic problems with the Izotov et al. (2007) determination related with 
the temperature structure are: (a) that to compute the A''(0++) abundance they adopted 
the T(4363/5007) value which is equivalent to adopt t'^ = 0.00, and (b) to compute the 
once ionized oxygen and nitrogen abundances, A^(0+) and A^(N+), they adopted the 
T(4363/5007) value, but according to photoionization models and observations of O- 
poor objects, the temperature in the 0+ regions, T(0+), is considerably smaller than 
in the 0++ regions, typically by about 2000 K for objects with T(0++) = 16000 K, 
and reaching 4000 K for the metal poorest H ii regions (e. g. Peimbert, Peimbert, & 
Luridiana 2002; Stasinska 1990). 

2.1. Recombination coefficients of the helium I lines 

There arc two recent line cmissivity estimates to derive the He abundance from recom- 
bination lines: one due to Benjamin, Skillman, & Smits (1999, 2002), and another by 
Bauman et al. (2005) and Porter et al. (2007, 2009). The difference in the Y values de- 
rived from both sets of data amoimt to about 0.0040, the emissivities by the first group 
yielding values smaller than those of the second group. According to Porter, Ferland, & 
MacAdam (2007), the error introduced in their emissivities by interpolating in temper- 
ature the equations provided by them is smaller than 0.03%, which translates into an 
error in Yp considerably smaller than 0.0001. Moreover according to Porter et al. (2009) 
the expected error in their line emissivities amounts to about 0.0010 in the Y derived 
values. In this review we are adopting the Bauman et al. and Porter et al. emissivities in 
the presented Yp values. 

2.2. Beyond case B 

There are at least four processes that modify the level populations of H and He atoms 
relative to case B and consequently the Yp determination: the optical depth of the He i 
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Table 1. COSMOLOGICAL PREDICTIONS BASED ON SBBN AND OBSERVATIONS 

FOR Tn = 885.7 ± 0.8 s 



Method 


Yp 


Dp 


?7io 




Yp 


0.2477 ± 0.0029" 


2.93 + 2.53 - LOe*" 


5.625 ± 1.81'' 


0.02054 ± 0.00661'' 


Dp 


0.2479 ± 0.0007' 


2.82 ± 0.28" 


5.764 ± 0.360'' 


0.02104 ± 0.00132'' 


WMAP 


0.2487 ± O.OOOe*" 


2.49 ± O.IS*" 


6.226 ± 0.170'' 


0.02273 ± 0.00062" 



"Observed value. ''Predicted value. References: t„ Arzumanov et al. (2000); Yp Peimbert et al. 
(2007a); Dp O'Meara etaL(2006); WMAP Dunkley et al. (2009). 



Table 2. COSMOLOGICAL PREDICTIONS BASED ON SBBN AND OBSERVATIONS 

FOR T„ = 878.5 ± 0.8 s 



Method 


Yp 


Dp 


Vio 




Yp 


0.2477 ± 0.0029" 


2.22 + 1.46 - 0.71'' 


6.688 ± 1.81'' 


0.02442 ± 0.00661'' 


Dp 


0.2462 ± 0.0007'' 


2.82 ± 0.28" 


5.764 ± 0.360'' 


0.02104 ± 0.00132'' 


WMAP 


0.2470 ± 0.0006'' 


2.49 ± 0.13'' 


6.226 ± 0.170'' 


0.02273 ± 0.00062" 



"Observed value. ''Predicted value. References: same as in Table 1 with the exception of Tn 
that comes from Serebrov et al. (2005, 2008). 

lines, the collisional excitation from the He 2 metastable level, the collisional excita- 
tions from the H ground level, and the fluorescent excitation of the H I and the He i 
lines, (case D of Luridiana et al. 2009). The last two processes are the least studied of 
the four. 

Luridiana et al. (2009) have introduced case D into the study of gaseous nebulae. Case 
D increases slightly the emissivities of the H i and He i lines affecting the accuracy of 
the Yp determination. There are no published estimates of the importance of this effect 
but it depends on the spectra of the ionizing stars, particularly in the region of the H I 
and He i lines, the spatial distribution of the ionizing stars, the gaseous electron density 
distribution, the fraction of ionizing photons that escape the nebula, the radial velocity of 
the gas relative to that of the stars, the nebular turbulence, and the region of the nebula 
observed. Therefore it requires tailor-made models for each observed object. Peimbert 
et al. (2007a) presented a list of thirteen sources of error in the Yp determination, each 
of them with systematic and statistical components but dominated by one or the other. 
The systematic error produced by not having considered Case D should be added to that 
list. 



3. Comparison of the directly determined Yp with the Yp values 
computed under the assumption of the SBBN and the 
observations of Dp and WMAP 

To compare the Yp value with the primordial deuterium abundance Dp (usually ex- 
pressed as 10^{D/H)p) and with the WMAP results, we will use the framework of the 
SBBN. The ratio of baryons to photons multiplied by 10^^,r]iQ, is given by (Steigman 
2006, 2007): 

7?io = (273.9 ± 0.3)nb/i^ (3.1) 

where Sib is the baryon closure parameter, and h is the Hubble parameter. In the range 
4 < r]io < 8 (corresponding to 0.2448 < 1^ < 0.2512), Yp is related to r]io by (Steigman 
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2006, 2007): 

Yp = 0.2483 ± 0.0005 + 0.0016(7?io - 6). (3.2) 

In the same rjio range, the primordial deuterium abundance is given by (Steigman 
2006, 2007): 

10^{D/H)p = Dp = 46.5(1 ± 0.03)(ryio)"^ ^ (3.3) 

From the Yp value by Peimbert et al. (2007a), the Dp value by O'Meara et al. (2006), 
the Ofc/i^ value by Dunkley et al. (2009), and the previous equations we have produced 
Table 1. From this table, it follows that within the errors Yp, Dp, and the WMAP 
observations are in very good agreement with the predicted SBBN values. 

Equations (2), (3), and (4) were derived under the assumption of a neutron lifetime, t„ 
, of 885.7 ±0.8 s (Arzumanov et al. 2000). A recent result by Serebrov et al. (2005, 2008) 
yielded a t„ of 878.5 ± 0.7 ± 0.3 s. This result would lead to a SBBN Yp value of 0.2470 
for the WMAP Ofe/i^ determined value (Steigman 2007, Mathews et al. 2005). Mathews 
et al. obtain for t„ = 881.9 ±1.6 s, mainly the average of the results by Arzumanov 
et al. (2000) and Serebrov et al. (2005, 2008), a Yp value 0.0009 smaller than for t„ = 
885.7 ±0.8. 

The 9(7 difference between both t„ determinations probably indicates that at least one 
of them includes systematic errors that have not yet been sorted out. 

From the Yp by Peimbert et al. (2007a), the Dp by O'Meara et al. (2006), and SBBN it 
is found that for t„ = 881.9 ± 1.6 the number of effective neutrino families, N^ff, is equal 
to 3.12 ± 0.23. Based on the production of the Z particle by electron-positron collisions 
in the laboratory and taking into account the partial heating of neutrinos produced by 
electron-positron annihilations during SBBN Mangano et al. (2002) find that N^ff = 
3.04. The Neff value derived from Yp, Dp, and SBBN is in excellent agreement with the 
value derived by Mangano et al. (2002). 

The Yp derived from WMAP is based on the very strong assumption of SBBN. It is 
also possible to derive Yp from the study of the microwave radiation without assuming 
SBBN: from the cosmic microwave background radiation {i.e. WMAP + ACBAR + 
CBI + BOOMERANG), Ichikawa, Sc^kiguchi and Takahashi (2008a,b) obtain that Yp < 
0.44, when they also include the information obtained from BAO + SN + HST (baryon 
acoustic oscillations in the distribution of galaxies, the distance measurements from type 
la supernovae, and the HST value for Hq) the constrain improves to Yp = 0.25 ^p ^y. 
Including the expected data from the Planck satellite they predict a reduction on the 
Yp error of about a factor of four to seven, an error still about four to six times higher 
than the one estimated from the best Yp determinations based on metal poor H II region 
observations. 

4. The AY/ AO ratio 

To determine the Yp value from a set of metal poor H ii regions it is necessary to 
estimate the fraction of helium, present in the interstellar medium of the galaxy where 
each H II region is located, produced by galactic chemical evolution. Prom observations of 
metal poor extragalactic H II regions it has been found that the Y versus O observations 
can be fitted with a straight line given by AY/ AO and equation (1) has been used often 
to derive Yp. A straight line is predicted by chemical evolution models of metal poor 
galaxies with the same initial mass function, a given set of stellar yields, and different 
star formation rates. To obtain different AY /AO from the models it is necessary to 
change the initial mass function (for example the maximum mass allowed or the slope at 
high masses), or the adopted yields. 
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The observational Yp, AF, and AO are affected by different amounts in the presence of 
temperature variations, while Yp diminishes by 0.0046 due to temperature variations in 
the sample of Pcimbcrt ct al. (2GG7a), Ay is slightly affected and AO is strongly affected 
by temperature variations (e. g. Carigi & Peimbert 2008, Table 3). In addition a fraction 
of O is embedded in dust grains, fraction that needs to be estimated to compare with 
models of galactic chemical evolution. 

The importance of Ay/ AO is two fold: it permits us to obtain a more accurate Yp 
value, and permits us to test for the presence of large temperature variations in gaseous 
nebulae when comparing nebular values with stellar ones. 

4.1. The gaseous 0/H determination 

There are two methods to derive the gaseous 0/H ratio, from the /(4363)//(5007) ratio 
together with the 7(3727) //(H/3) and the /(5007)/7(H/3) line ratios, the so called T(4363) 
method, and from the intensity ratio of O II recombination lines to H i recombination 
lines that has been called the O ii^l method by Peimbert et al. (2007b). The O Url 
method usually provides higher 0/H ratios by 0.15 to about 0.3 dex, this difference is 
due to temperature variations inside the observed volume and is smaller for metal poor 
H II regions and higher for metal rich H ii regions. The T(4363) method in the presence 
of temperature variations produces a systematic effect that lowers the 0/H abundances 
relative to the real ones. On the other hand the O Url method is independent of the 
temperature structure. 

4.2. The total 0/H determination 

Another factor that has to be taken into account to obtain the total 0/H ratio is the 
fraction of O atoms trapped in dust grains. Esteban et al. (1998), based on the depletion 
of Fe, Mg and Si in the Orion nebula, estimated that the fraction of oxygen atoms trapped 
in dust grains amomits to 0.08 dex. Mesa-Delgado et al. (2009) have estimated that the 
fraction of O atoms trapped in dust grains in the Orion nebula amounts to 0.12 ± 0.03 
dex; this result is based on three different methods: a) by comparing the O abundances of 
the B stars of the Orion association with the O abundance of the Orion nebula, b) from 
the depletion of Fe, Mg and Si in the nebula and assuming that these atoms are combined 
with molecules containing O, and c) by comparing the shock and nebular abundances of 
the material associated with the Herbig-Haro object 202. 

Rodriguez & Rubin (2005) have estimated the fraction of Fe atoms in galactic and 
extragalactic H ii regions in the gaseous phase, the depletions derived for the different 
objects define a trend of increasing depletion at higher metallicities; while the Galactic 
H II regions show less than than 5% of their Fe atoms in the gas phase, the extragalactic 
ones (LMC 30 Doradus, SMC N88A, and SBS 0335-052) have somewhat lower depletions. 
Izotov et al. (2006) find a slight increase of Ne/0 with increasing mctallicity, which they 
interpret as due to a moderate depletion of O onto grains in the most mctal-rich galaxies, 
they conclude that this 0/Ne depletion corresponds to ~ 20% of oxygen locked in the 
dust grains in the highcst-mctallicity H II regions of their sample, while no significant 
depletion would be present in the H ii regions with lower metallicity. Peimbert & Peimbert 
(2010) based on the Fe/0 abundances of Galactic and extragalactic H ii regions estimate 
that for objects in the 8.3 < 12 + log 0/H < 8.9 range the fraction of O atoms trapped 
by dust grains amounts to 0.12 ± 0.03 dex, for objects in the 7.7 < 12 + log 0/H < 8.3 
range amounts to 0.09±0.03 dex, and for objects in the 7.3 < 12 + log 0/H < 7.7 range 
amounts to 0.06 ± 0.03 dex. 
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4.3. Extrapolation of the Y determinations to the value ofYp, or the O (AY/ AO) 

correction 

From chemical evolution models of different galaxies it is found that AY/ AO depends 
on the initial mass function (IMF), the star formation rate, the age, and the O value of 
the galaxy in question. Peimbert et al. (2007b) have found that AY/ AO is well fitted 
by a constant value for objects with the same IMF, the same age, and an O abundance 
smaller than 4 x 10~^. This result is consistent with the custom of using a constant value 
for AY/ AO to fit observational data. 

To obtain an accurate Yp value, a reliable determination of AY/ AO for 0-poor ob- 
jects is needed. The AY/ AO value derived by Peimbert, Peimbert, & Ruiz (2000) from 
observational results and models of chemical evolution of galaxies amounts to 3.5 ± 0.9. 
More recent results are those by Peimbert (2003) who finds 2.93 ±0.85 from observations 
of 30 Dor and NGC 346, and by Izotov & Thuan (2004) who, from the observations of 
82 H II regions, find AY/ AO = 4.3 ± 0.7. Peimbert et al. (2007a) have recomputed this 
value by taking into account two systematic effects not considered by Izotov & Thuan: 
the fraction of oxygen trapped in dust grains, and the increase in the O abundances due 
to the presence of temperature variations. Prom these considerations they obtained for 
the Izotov & Thuan sample that AY/ AO = 3.2 ± 0.7. 

On the other hand Peimbert et al. (2007b) from chemical evolution models with dif- 
ferent histories of galactic inflows and outflows for objects with O < 4 x 10~^ find that 
2.4 < AY/ AO < 4.0. From the theoretical and observational results Peimbert et al. 
(2007a) adopted a value of AY/ AO = 3.3 ± 0.7, that they used with the Y and O 
determinations from each object to obtain the Yp value. 

4.4. Comparison of the oxygen abundances of the ISM of the solar vicinity with those of 
the Sun and F and G stars of the solar vicinity 

In addition to the evidence presented in section 2 in favor of large t"^ values, and con- 
sequently in favor of the O Url method, there is another independent test that can be 
used to discriminate between the r(4363) method and the O Ubl method that consists 
in the comparison of stellar and H II region abundances of the; solar vicinity. 

Esteban et al. (2005) determined that the gaseous 0/H value derived from H II regions 
of the solar vicinity amounts to 12 -I- log (0/H) = 8.69, and including the fraction of O 
atoms tied up in dust grains it is obtained that 12 -I- log (O/H) = 8.81 ± 0.04 for the 
O/H value of the ISM of the solar vicinity. Alternatively from the protosolar value by 
Asplund et al. (2009), that amounts to 12 + log(0/II) = 8.71, and taking into account the 
increase of the O/H ratio due to galactic chemical evolution since the Sun was formed, 
that according to the chemical evolution model of the Galaxy by Carigi et al. (2005) 
amounts to 0.13 dex, we obtain an 0/H value of 8.84 ± 0.04 dex, in excellent agreement 
with the value based on the O Wrl method. In this comparison we are assuming that 
the solar abundances are representative of the abundances of the solar vicinity ISM when 
the Sun was formed. 

There are two other estimates of the 0/H value in the ISM that can be made from 
observations of F and G stars of the solar vicinity. According to Allende-Prieto et al. 
(2004) the Sun appears deficient by roughly 0.1 dex in O, Si, Ca, Sc, Ti, Y, Ce, Nd, and 
Eu, compared with its immediate neighbors with similar iron abundances; by assuming 
that the O abundances of the solar immediate neighbors are more representative of the 
local ISM than the solar one, and by adding this 0.1 dex differencx; to the protosolar solar 
value by Asplund et al. (2009) we obtain that the present value of the ISM has to be 
higher than 12 -|- log 0/H = 8.81. A similar result is obtained from the data by Bensby 
& Feltzing (2005) who obtain for the six most 0-rich thin-disk F and G dwarfs of the 
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solar vicinity an average [0/H] — 0.16; by assuming their value as representative of the 
present day ISM of the solar vicinity we find 12 + log 0/H = 8.87. Both results are in 
good agreement with the 0/H value derived from the O Url method. 

4.5. Comparison of /SY//S.Z values of the Galactic Hu region M17 with K dwarfs of 

the solar vicinity 

The best Galactic H ii region to determine the He/H ratio is M17 because it contains 
a very small fraction of neutral helium and the error introduced by correcting for its 
presence is very small. Carigi & Peimbert (2008) obtained from observations for M17 
(and adopting a Yp value) a value of AY/AZ = 1.97 ± 0.41 for = 0.036 ± 0.013, 
where Y and Z are the helium and heavy elements by unit mass; by correcting this value 
considering that the fraction of O trapped in dust amounts to 0.12 dex instead of 0.08 
dex (Mesa-Delgado et al. 2009) we obtain AY/AZ = 1.77 ± 0.37. This AY/AZ value is 
in agreement with two independent AY/AZ determinations derived from K dwarf stars 
of the solar vicinity that amount to 2.1 ± 0.4 (Jimenez et al. 2003) and to 2.1 ± 0.9 
(Casagrande et al. 2007). 

On the other hand the value AY/AZ = 3.6 ± 0.68 derived from coUisionally excited 
lines of M17 under the assumption of t'^ = 0.00 by Carigi and Peimbert (2008) (value 
corrected for a fraction of O trapped by dust grains of 0.12 dex) is not in agreement with 
the values derived from K dwarf stars of the solar neighborhood. 

4.6. Comparison of Yp and AY/ AO with models of chemical evolution for the disk of 

the Calaxy 

To compare Galactic chemical evolution models of Y and O with observations we need to 
use the best determinations available of the abundances of these elements. We consider 




Figure 1. In the two panels the origin corresponds to 1^ = 0.2477, the triangles represent the 
Y and O values derived by Peimbert et al. (2007a) for five metal poor extragalactic H ll regions, 
the open and filled circles represent the M17 values for t^ = 0.00, and = 0.036 respectively 
from Carigi & Peimbert (2008) , the large open circle with a dot in the center corresponds to the 
presolar values by Asplund et al. (2009). The left panel presents two chemical evolution models 
by Carigi & Peimbert (2010), the solid line is a model for the Galactic disk with a time span 
from the formation of the Galaxy to the formation of the Sun at a galactocentric distance of 8 
kpc, while the dashed line corresponds to a model for the Galactic disk with a time span from 
the formation of the Galaxy to the present at a galactocentric distance of 6.75 kpc, the distance 
of M17 to the galactic center; the stellar yields for the two panels are different, see subsection 
4.6 for further details. 
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that the two most accurate Galactic Y and O determinations are the presolar values 
(Asplund et al. 2009) and the M17 H ii region values (Carigi & Peimbert 2008). 

Carigi et al. (2005) presented chemical evolution models for the disk of the Galaxy that 
fit the slope and the absolute value of the 0/H gradient, these models are also successful 
in reproducing the C/0 gradient derived from H II regions and the C/0 versus 0/H 
evolution history of the solar vicinity obtained from stellar observations. In Figure 1 
we present chemical evolution models by Carigi & Peimbert (2010) for the Y and O 
abundances of the Galactic disk for two sets of stellar yields where they have added the 
presolar values by Asplund et al. (2009) for comparison. These models are based on the 
same assumptions than those adopted by Carigi and Peimbert (2008). The only difference 
between the models in the left panel and those in the right one, is that for massive stars 
with Z > 0.004 the ones in the left use the yields by Maeder (1992) with high mass loss, 
while the ones on the right use the yields by Hirschi, Meynet, & Maeder (2005) with 
low mass loss. There are at least three conclusions that we can extract from the figure: 
(a) the fit for the presolar values is very good for the two sets of yields, (b) the fit for 
M17 for = 0.036 ± 0.013 is good, but it can be improved by assuming a set of yields 
intermediate between the two sets used; a similar result in favor of an intermediate set 
of yields was obtained by Cescutti et al. (2009) based on C/0 observations for stars in 
the bulge of the Galaxy, (c) the models do not fit M17 for t'^ = 0.00, this result is in 
agreement with those by Esteban et al. (2005, 2009) and Peimbert et al. (2007a) who 
find large values for Galactic and extragalactic H ii regions. Moreover the O and Y 
abundances for the presolar material and for M17 arc derived from independent methods, 
and they are fitted by the same chemical evolution model. This fit provides us with a 
consistency check on the gaseous nebulae helium and oxygen abundance determinations 
based on large t^ values. 

5. Conclusions 

During the last 50 years the determination of Yp has been very important for the study 
of cosmology, stellar evolution, and the chemical evolution of galaxies. To determine^ Yp it 
is necessary to determine accurate atomic parameters and the physical conditions inside 
ionized gaseous nebulae. 

During the last five decades the accuracy of the Yp determination has increased consid- 
erably, and during the last two decades the differences among the best Yp determinations 
have been due to systematic effects (e. g. Olive & Skillman 2004, Peimbert 2008). These 
systematic effects have been gradually understood, particularly during the last few years. 

The best Yp determination available, that by Peimbert et al. (2007a), is in agreement 
with the Dp determination and with the WMAP observations under the assumption of 
SBBN. The errors in the Yp determination are still large and there is room for non- 
standard physics. 

Similarly to improve the accuracy of the Yp value derived under the assumption of 
SBBN and the Qbh"^ derived from the background radiation, a new determination of 
the neutron lifetime is needed to sort out the difference between the r„ obtained by 
Arzumanov et al. (2000) and the t„ obtained by Serebrov et al. (2005, 2008). 

To improve the accuracy of the Yp determination based on metal poor H II regions, the 
following steps should be taken in the near future: (a) to obtain new observations of high 
spectral resolution of metal poor H II regions, those with 0.0005 < Z < 0.001 to reduce 
the effect of the coUisional excitation of the Balmer lines, which for the present day Yp 
determinations is one of the two main sources of error; (b) to determine the temperature 
of the H II regions based on a large number of He i lines observed with high accuracy, 
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and to try to avoid the use of T(4363/5007) temperatures that weigh preferentiaUy the 
regions of higher temperature than the average one, effect that artificially increases the 
Y.p determinations; (c) additional efforts should be made to understand the mechanisms 
that produce temperature variations in giant H ii regions and once they are understood 
they should be incorporated into photoionization models; (d) the He I recombination 
coefficients should be computed again with an accuracy higher than that of the last two 
determinations; (e) the computation of tailor-made photoionization models for extra- 
galactic H II regions including the effect of case D on the intensity of the H i and He i 
lines should be carried out. 

The observed AY/ AO ratios provide us with strong constrains for models of galactic 
chemical evolution. The gaseous O abundances have to be corrected by the fraction of 
oxygen embedded in dust grains and by the effect of temperature variations, these two 
corrections together amount to about 0.2 to 0.4 dex. Chemical evolution models for the 
Galactic disk are able to reproduce the observed Y and O presolar values and the Y and 
O values derived for M17 based on H, He and O recombination lines, but not the M17 
Y and O values derived from r(4363/5007) and O coUisionally excited lines under the 
assumption of t'^ = 0.00. This result provides a consistency check in favor of the presence 
of large temperature variations in H II regions and in favor of the Yp determinations 
based on the T(He i) and values derived from He I recombination lines. 
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